{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "8ca704e1",
   "metadata": {},
   "source": [
    "# 矩形窗"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "dee9b079",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAESCAYAAAAPEjVtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAASJUlEQVR4nO3dfZBdd13H8feneXDWVGhql9h2hFqt1Uop1IUSLXXp9IkZeTAKjWIdniYyVvFpqkWCOk6hTGU6CjNliJaKikBBGlSoBB8CVduBzYTypI6ILRjIsNiH2MoIlq9/3Jtmu7l378nmbrL72/drJpNzz/nec7+/Xzafe+acc++mqpAkrWwnHO8GJElHzzCXpAYY5pLUAMNckhpgmEtSAwxzaZGSrE+y5nj3IYFhrsYk2ZjkXUm+pUPt2gHr1s9ZflqSb19gF7uACxfXqTReh/0wSyvcJcDXq+p/oRfuwOeBf59TswH4FeBVSS4Avkzv/8LJwOeAzUkeB9wG/DJwWz/k/2neaz0R+OMks/PW/0hVPTzeYUkLM8zVmhcDFyT5PL2f72uAh6tqakDt7UneBbwBOAnYWlWvTHICcAvwB1V1W7/2EeB7quqkJD8GPK2qfrN/dP8e4MPAO4CPAFnC8UkDeZpFzUhyLvC9wI8AXwC+G3jfEe5jDfBH9N4AXpdkXX9T0TsK/zhwKfD+JLcD6+gd5W8C9tJ7E/CoXMecR+ZqyQHgZcA3AarqGwBJnpDkE/2aE4FPAVvoBfF8JwDvBj6UZCvw0iRXVNU3k9wC/A/wVeCXgI3A1f3n/Te9N4EPl9+RoePAMFdLzgY2A48HzkxyA/A3wFeq6qkASS4BrgK+C/gAcDrwFGA9vVMtzwKeDkwD1wGXzQnnLwJXAi+l98YxC8zMef3fo3fKRTrmDHO1ZBb4OPB9wD7g7+hd3DxMVX0e+P4k/0jvKP37ga3AzwE/D/wicEW/7qD/o3cK5bPAm4BvB87vb3sb8ADwjbGOSOrIc+ZqRlXtraoP0vu53gTsrapPAU9IMpNkht7Fzn1znnYqvTeBg54E/FD/z/4kf5vkfB7r9f393AJso3cXzAeWYkxSVx6Zqyn9O1FeQu9i5IeTXETvNMthd7MkeTLwpf75cACq6j+ArUlOBXb3/3xi3lN/GTgP+A3g14ELq+q/Du5DOh48MldrrgceBF4I3AncML8gybok3wr8GvD2/uoT5mx/PnAXcFNV/WpVfbO/6eDBzzTwqv5rXUvvKJ75+5GOpXjhXa1IciVwI3BBVf1nkgngCfTu/57rBOBD9G4xvLiqvp7k6fRuaZwBfh94SVXtnbf/U+l9AOkzQ1o4G9hcVZ8e15ikrgxzNSO98xynVtWXjnY/3l6olcYwl6QGeH5PkhpgmEtSA47LrYmnnHJKnXHGGcfjpSVpxdqzZ89Xq2py0LbjEuZnnHEGMzMzowslSY9Kcu+wbZ5mkaQGGOaS1ADDXJIaYJhLUgMMc0lqgGEuSQ0wzCWpAYa5JDXAMJekBhjmktSATmGeZFOSOzrW7R1VJ0kar5FhnmQjvV+ttaHD/t4ITBxtU5KkI9PlyPwR4ErgwEJFSS4GHgb2j6EvSdIRGBnmVXWgqh5cqCbJeuC19H657bCabUlmkszMzs4eeaeSpKHGdQH0Wnq/yfyBYQVVtaOqpqpqanJy4NfxSpIWaVxhfglwdZLdwFOT/OGY9itJ6uCIfzlFknOAn6qq7QfXVdVFc7bvrqpXjKk/SVIHncO8qqb7f38W2D6qTpJ07PihIUlqgGEuSQ0wzCWpAYa5JDXAMJekBhjmktQAw1ySGmCYS1IDDHNJaoBhLkkNMMwlqQGGuSQ1wDCXpAYY5pLUAMNckhpgmEtSAwxzSWqAYS5JDTDMJakBhrkkNcAwl6QGdArzJJuS3LHA9scnuT3JriS3JVk/vhYlSaOMDPMkG4G3AxsWKHsxcGNVXQbsB64YT3uSpC7Wdqh5BLgSeP+wgqq6ac7DSeArR9mXJOkIjDwyr6oDVfVgl50l2QxsrKq7BmzblmQmyczs7OwiWpUkDTO2C6BJTgbeDLxs0Paq2lFVU1U1NTk5Oa6XlSQxpjDvX/B8D/Dqqrp3HPuUJHV3xGGe5Jwk181b/XLgfOA1SXYnuXIs3UmSOklVHfMXnZqaqpmZmWP+upK0kiXZU1VTg7b5oSFJaoBhLkkNMMwlqQGGuSQ1wDCXpAYY5pLUAMNckhpgmEtSAwxzSWqAYS5JDTDMJakBhrkkNcAwl6QGGOaS1ADDXJIaYJhLUgMMc0lqgGEuSQ0wzCWpAYa5JDXAMJekBqztUpRkE/DeqnrWAjU3A+cAH6iq68bU36N27t3H737oX/nSA1/jtJMmuObys3nB005fcJvrnSPnwjla7nMxLqmqhQuSjcA7gSdU1flDarYAz6uqlyR5G3B9Vf3bsH1OTU3VzMxM5yZ37t3Hq9/3Kb72jUceXTexbg3XbzkXYOC2H//B0/nzPftW/XrnyLlwjpbvXFy/5dwjCvQke6pqauC2DmH+OCDA+6tqekjNm4C/rqoPJtkKTFTVLcP2eaRh/sNv+Dv2PfC1w9afftIEwMBtaxIeGTC21bbeOTrEuRi93jk65FjMxeknTfCP11582PphFgrzkadZqupAfycLlW0A9vWX7wMOO4JPsg3YBvDEJz5x1Ms+xpcGTOZC64GBE7ca1ztHhzgXo9c7R4cci7lY6DWO1LgugD4ETPSXTxy036raUVVTVTU1OTl5RDs/7aSJoeuHbVsz5M1nta13jg5xLkavd44OORZzMWz/izGuMN8DXNhfPg+4Z0z7BeCay89mYt2ax6ybWLeGay4/e+i2n7zgO13vHDkXztGynotrLj+bcel0N8tcSc4Bfqqqts9ZvRO4I8lpwHOAZ46nvZ6DFwh+7b2f5OuPfJPTB1wJHrRt6kknu945ci6co2U9F+My8gJo5x317nq5FPhoVe1fqPZIL4AedOVb7wTg3T+7ufM21x//115u65djT8tt/XLsqeW56OqoLoB2VVX3A7eOa3+SpO78BKgkNcAwl6QGGOaS1ADDXJIaYJhLUgMMc0lqgGEuSQ0wzCWpAYa5JDXAMJekBhjmktQAw1ySGmCYS1IDDHNJaoBhLkkNMMwlqQGGuSQ1wDCXpAYY5pLUAMNckhrQKcyT3JzkziTbh2zfmOSDSWaSvHW8LUqSRhkZ5km2AGuqajNwZpKzBpRdBbyjqqaAb0syNeY+JUkL6HJkPg3c2l/eBVw4oOa/gCcnOQn4TuCL42hOktRNlzDfAOzrL98HbBpQ8w/Ak4BXAf/cr3uMJNv6p2FmZmdnF9muJGmQLmH+EDDRXz5xyHN+C3hlVf0O8C/AS+cXVNWOqpqqqqnJycnF9itJGqBLmO/h0KmV84B7BtRsBM5Nsga4AKixdCdJ6qRLmO8ErkpyI/Ai4DNJrptXcz2wA3gQOBl45ziblCQtbO2ogqo6kGQauBS4oar2A3fPq/kY8ANL0aAkabSRYQ5QVfdz6I4WSdIy4ydAJakBhrkkNcAwl6QGGOaS1ADDXJIaYJhLUgMMc0lqgGEuSQ0wzCWpAYa5JDXAMJekBhjmktQAw1ySGmCYS1IDDHNJaoBhLkkNMMwlqQGGuSQ1wDCXpAYY5pLUgE5hnuTmJHcm2T6i7qYkzx1Pa5KkrkaGeZItwJqq2gycmeSsIXXPAr6jqv5yzD1KkkbocmQ+DdzaX94FXDi/IMk64A+Ae5I8f2zdSZI66RLmG4B9/eX7gE0Dan4G+CxwA/CMJL8wvyDJtiQzSWZmZ2cX268kaYAuYf4QMNFfPnHIc54G7Kiq/cCfAs+eX1BVO6pqqqqmJicnF9uvJGmALmG+h0OnVs4D7hlQ8zngzP7yFHDvUXcmSepsbYeancAdSU4DngNsTXJdVc29s+Vm4G1JtgLrgJ8Ye6eSpKFGhnlVHUgyDVwK3NA/lXL3vJr/Bl64FA1KkkbrcmROVd3PoTtaJEnLjJ8AlaQGGOaS1ADDXJIaYJhLUgMMc0lqgGEuSQ0wzCWpAYa5JDXAMJekBhjmktQAw1ySGmCYS1IDDHNJaoBhLkkNMMwlqQGGuSQ1wDCXpAYY5pLUAMNckhpgmEtSAzqFeZKbk9yZZPuIuk1J9o6nNUlSVyPDPMkWYE1VbQbOTHLWAuVvBCbG1ZwkqZsuR+bTwK395V3AhYOKklwMPAzsH0tnkqTOuoT5BmBff/k+YNP8giTrgdcC1w7bSZJtSWaSzMzOzi6mV0nSEF3C/CEOnTo5cchzrgVuqqoHhu2kqnZU1VRVTU1OTh5xo5Kk4bqE+R4OnVo5D7hnQM0lwNVJdgNPTfKHY+lOktTJ2g41O4E7kpwGPAfYmuS6qnr0zpaquujgcpLdVfWKsXcqSRpqZJhX1YEk08ClwA1VtR+4e4H66XE1J0nqpsuROVV1P4fuaJEkLTN+AlSSGmCYS1IDDHNJaoBhLkkNMMwlqQGGuSQ1wDCXpAYY5pLUAMNckhpgmEtSAwxzSWqAYS5JDTDMJakBhrkkNcAwl6QGGOaS1ADDXJIaYJhLUgMMc0lqgGEuSQ3oFOZJbk5yZ5LtQ7Y/PsntSXYluS3J+vG2KUlayMgwT7IFWFNVm4Ezk5w1oOzFwI1VdRmwH7hivG1KkhaytkPNNHBrf3kXcCHwb3MLquqmOQ8nga+MozlJUjddTrNsAPb1l+8DNg0rTLIZ2FhVdw3Yti3JTJKZ2dnZRTUrSRqsS5g/BEz0l08c9pwkJwNvBl42aHtV7aiqqaqampycXEyvkqQhuoT5HnqnVgDOA+6ZX9C/4Pke4NVVde/YupMkddIlzHcCVyW5EXgR8Jkk182reTlwPvCaJLuTXDneNiVJCxl5AbSqDiSZBi4Fbqiq/cDd82reArxlKRqUJI3W5W4Wqup+Dt3RIklaZvwEqCQ1wDCXpAYY5pLUAMNckhpgmEtSAwxzSWqAYS5JDTDMJakBhrkkNcAwl6QGGOaS1ADDXJIaYJhLUgMMc0lqgGEuSQ0wzCWpAYa5JDXAMJekBhjmktQAw1ySGtApzJPcnOTOJNuPpkaStDRGhnmSLcCaqtoMnJnkrMXUSJKWTqpq4YLkTcBfV9UHk2wFJqrqliOtmWtqaqpmZmaOuNlnvv5vePjrj3DOqY87bNtnv3wA4LBtrj/+r73c1i/Hnpbb+uXYUytzsWH9Gu76jUsOe40ukuypqqmB2zqE+c3Am6rq7iSXAedX1RsWUbMN2NZ/eDbwr4saDZwCfHWRz12pHPPq4JhXh6MZ85OqanLQhrUdnvwQMNFfPpHBp2ZG1lTVDmBHh9dbUJKZYe9MrXLMq4NjXh2WasxdLoDuAS7sL58H3LPIGknSEulyZL4TuCPJacBzgK1Jrquq7QvUPHPcjUqShht5ZF5VB4Bp4C7g2VV197wgH1Tz4PhbfdRRn6pZgRzz6uCYV4clGfPIC6CSpOXPT4BKUgMMc+k4SXJykkuTnHK8e9HKt6LCfDV9ZUCSTUnumPO46bEneXyS25PsSnJbkvUtjznJRuCvgGcAf59ksuXxztX/2d7bX256zEnWJvlCkt39P+cu1ZhXTJivpq8M6P9Hfzuwof94NYz9xcCNVXUZsB/YSttjfgrwK1X1OuBDwMW0Pd653ghMrJKf66cA76yq6aqaBs5iica8YsKc3t0yt/aXd3HovvYWPQJcCRzoP56m8bFX1U1V9eH+w0ngp2l4zFX1kaq6K8lF9I7OL6fh8R6U5GLgYXpv2NO0P+ZnAj+a5GP9T8pfwhKNeSWF+QZgX3/5PmDTcexlSVXVgXm3d66asSfZDGwEvkjjY04Sem/a9wNF++NdD7wWuLa/ajX8XH8cuKSqngGso/c5nCUZ80oK8y5fK9CqVTH2JCcDbwZexioYc/VcDXwS+CEaHy+9EL+pqh7oP27+3xj4ZFV9ub88Q+97WZZkzCtp8lbzVwY0P/b+Udt7gFdX1b00PuYkv57kZ/oPTwLeQMPj7bsEuDrJbuCpwHNpf8x/kuS8JGuAFwBXs0Rj7vJx/uViJ6v3KwN20v7YXw6cD7wmyWuAW4CrGh7zDuDWJK8APk3v3/ijDY+Xqrro4HI/0J9H+z/XvwP8GRDgL1jC/8sr6hOg/bs8LgU+WlX7j3c/x9JqHPtqG/NqGy845nGOeUWFuSRpsJV0zlySNIRhLkkNMMwlqQGGuSQ1YCXdmigtqSS/Te9Tes8CHgdcsVrusNDK55G59Fjf078f+n30vvxKWhEMc+mx/rj/9xeA9cezEelIGObSYz18vBuQFsMwl6QG+AlQSWqAR+aS1ADDXJIaYJhLUgMMc0lqgGEuSQ0wzCWpAYa5JDXg/wFs9nrKGP9IDgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEXCAYAAABsyHmSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABp50lEQVR4nO29d3hkV3n4/3mnqdeVtnfvrtZre9dlXdm1tcamhBLi0EJiAik2iUMa4QchhORLgBAgkEIgNtimY4wJzd3YK3vX9u56e++9aNXLSJp+fn/ce2dG0pR7r6SdkXU+z6NHM1f3zH3nzOi8961HlFJoNBqNRuMGT6EF0Gg0Gs3kRSsRjUaj0bhGKxGNRqPRuEYrEY1Go9G4RisRjUaj0bhGKxGNRqPRuEYrEY0mByISEBFvnnMkwzHfBMpUOVGvrdE4RSsRzaRDROpE5BERKbFx7qjFXEQCaY+vEZFpOV7iWWBNnsvcJCLPjTj2oohcm0Wm34jIipGypP3dn01xicibgKfzyDNyzAdE5AM2z10gIvc6eX3N1GbC7pY0mgnkDiCilAqDoVSA48CxtHMqgL8F/lJEbgQuYHzf64GjwM0iUg38HPgb4Ofmgv7KiGvNB74nIu0jjt+mlBowH18FvGz9QURqgQZgRxb5Q0DEfPyfInLViL8HgI8BG0Tk7cBn0s73AleKyEbzeQnwnFLqU6ZSfYd5vEsp9YL5+C7gR1lkGUkE+HcR+YlSqsdUZn6lVMjmeM0UQ3TFumayISK/AG4EhjAUw8eBf1dKzc1y/iPAF4Fa4P1KqY+IiAf4KbBdKfV58zwv0KmUqhWR3wGuUUp9xrRmfgo8B/wQeBFYo5QKisivgetNWXqBzwHVwL8C/aYIM5VSlaYy+y9gJobSuUwpdaOIPI+hHCy+rpR6zJTpA8BqpdTfikg9xoJ+0ZT/WqXU1rT3OQ14Ffgq8EHgEeDvgNMjpmS1UqrUHLMQ2AXsyTLdPkAppW7O8nfNFEdbIppJhXnXvgy4DXgAuNP80787eA0v8DAwoJT6vIj4lVJRQGFYHa8BrwEPishTGHfyfwt8CMO6+E/AskIuA5abd+2fA6qAPwRuVUodMq+XbiHtAOZiWBffNo/5lVK3mue+H2hMFxdQppXxHmCOORbgcRGZBfhM+SPAeaXU/5qvEwEeUkr984j33zdiSnZZ1zf/bs2HRpMXrUQ0k40+4I+ABIC12InIdBHZaZ5TiXFnfRfgz/AaHuAnwDPmYvthEXmLUiohIg8Dg0AH8NdAHXCfOa4f+A6G+8gy4RMjXvs2QyxDgYw4Z+S5FstE5Afm4wUYVo+FADGgBZgGBETEWvDrMKyijcCnMJRgOjEAEflX4M0YbrI3Z5Aj+VxE7sZQjP9sPn8I+D+l1ONZZNdMcbQS0Uw2moCbgRpgsYh8CfgN0KaUuhpARO4A7gYWAU9g3L2vxFhEa4G1GC6oZgz305vSlMIZ4H3AhzEUVjuQdBkB/8HwRX4kW4Bvi8iDSqk/No9Zi3Q23/HbMJSFxZm0x6VASCl1s4jcA9wEPGT+7dF0CyIHs4A/x3BtjVSqHiAiIl/DiA/9HNhpPg9ixFg+a+MamimKViKayUY7hqtpOXAOeAEjaD4KpdRx4HIReRnDKrkceD/GgvoXwF8BbzHPs4hhuKr2Y8QvpgFWltVDQA8w0tWzQUTiGLGOTyqlNorIV0TkzUqpZ9LOy5QK/E0MpRbBCOIfBFaKyNWmXBVAmxn034mRGGBli30xyxw5oQ7j/f4U+DqwGsNK+ytgPUbM6OQ4XEfzOkUrEc2kQim1A9hhpsjOAHaYgebpIpJuMTyb9ngWhvK53Hy+ALjF/Bk0A9sfV0ptTxvzBYxFegXwPeAxDKvmjzKItTYtJmLx78AHReRFjDt6GK1EPBgK6T+Ai6QU2yOksrHmYATLv0Qq88qyVG4UkXal1I/TXnOFmUgwMn3/fvO1RjIXOKWUekVEzgC/A3wLQ3HejRE70miyoutENJMOMzPpQxhB6ufMlNo2pdTqtJ9PmedeiRFsTvr9lVInlFLvx/j+t2Dc4e8ccZm/wchqehewD/hTpVSnAzF/jqFwajDcYjDclfTXwGEMV9dajASBOcC7Mawfi1XAUaXUX2NYIF3AezGC8r8YoUACwH7zvSVlVUp9SCm1SinVkMGqWGnKAYbiXGSe86wp088dvGfNFEQrEc1k5F8x0mnfQ+oufRhmwV458P8B3zUPe9L+/tvAJuAbSqmPpSkZyzpvBv7SvNYnMawWRr4OhmLYYAb1/4SUokhgWB5vwahLAdiMoZx8wH8ppf7EfK3rMBTETPP8GlPGGgzraR+AUuoC8E8YMZqPMdoqCgEPmuf+zgg5rer7hQyPzbwF2C0iAaXUZqXUV8zjvwf8TCkVNseNcsVpNKDdWZpJhoi8D/gD4EYzm+qvgekYbpyNaad6gGeAxRiLOxjZVUdFpBkjWPwu0z2Wjh8jO+mfzOd/k3bt38cI7KdXmX8Z+L652K7CqBcBw9X2PxiK6mMASqkYEBMRf9prBIC/U0olq9BFZJN5/K3Ar4A1ZtHhFRjFhV/DKHDcKiIXMVxhf6eUOgd8P022kRX9DRgxpK+b17kBw7VXDpwRkVFpvWb2WgmGRXR+5N81Gl1sqJlUmHfEs5RSY1rQRERUEXz5zar5kFIqkuXvfozFfxWwTSnVPuLvi4GlIwL4dq9djqGM1zuXXKMx0EpEo9FoNK7RMRGNRqPRuGbSxERE5EGMdMsnlFKfy3ZebW2tWrJkyaUTLAsDAwNUVFRMeRmKRY5ikKFY5CgGGYpFjmKQoVjk2LZtW4dSqjH/mSNQShX9D0ah2HfMxw9h+IAznrvgsmVqKBJThWb9+vWFFqEoZFCqOOQoBhmUKg45ikEGpYpDjmKQQanikAPYqlysz5PFndUMPGo+zrm/w4WBBNd/7jc8tPEEiYSO92g0Gs1EMikC66Yr67+UUrvMTXmuVUp9Me3v9wD3ANROn33dmo8/yJ6OOLfP93H35QEKkeIeDAaprCzsBnTFIEOxyFEMMhSLHMUgQ7HIUQwyFIsc69at26aUWu14oBvz5VL/YLTevkmlXFufynbusmXLVCKRUF94Yr9a8InH1XdePjFWK88VxWCeFoMMShWHHMUgg1LFIUcxyKBUcchRDDIoVRxy8Dp3Z20j5cJaBZzMdbKI8Mm3Lue2ZY184ckDnOocyHW6RqPRaFwyWZTIL4C7ReSrGH2Dnsg3QET40rtX4hHhK88ezne6RqPRaFwwKZSIUqoPI7i+CVinlOq1M25GdSl/vGYRv951noOtIzdz02g0Gs1YmRRKBEAp1a2UelQp1epk3J+sXUSJz8N3Xzk5QZJpNBrN1GXSKBG31JYHeNfVc/j5jnP0DuptozUajWY8ed0rEYDfv2k+oWiCp/dl3ABPo9FoNC6ZEkrkqjk1zK8v5/HdWoloNBrNeDIllIiI8PaVs3jlWCddAxk7bms0Go3GBVNCiQC8+YqZxBOKjUc7Ci2KRqPRvG6YMkrkyjk11JT52XC4Pf/JGo1Go7HFlFEiXo+wZkkDG450WK1UNBqNRjNGpowSAViztIHWvhDHO3QbFI1GoxkPppQSuW5BHQA7T/cUVhCNRqN5nTCllMhljZVUlvjYcaa70KJoNBrN64IppUS8HmHl3Bp2nukptCgajUbzumBKKRGAa+bXcuBCP0OReKFF0Wg0mknPlFMiV86uIZ5QHGnrL7QoGo1GM+kpKiUiIjNEZMOIYw+KyKsi8unxuMbSGVUAHLkYHI+X02g0milN0SgREakDvgtUpB27C/AqpW4GFovI0rFeZ+G0cvxe4bC2RDQajWbMSLEU3olINSDAL5VSzeax/wKeVko9KSLvB8qUUg9nGHsPcA9AY2PjdY8++mjOa3164yDTyjz8zXWl4/wuUgSDQSorKyfs9SeLDMUiRzHIUCxyFIMMxSJHMchQLHKsW7dum1JqteOBbjZmH48f4H6gJe3nM+bxlrRzHgRWmY/fBHwy3+suW7Ys74b09/1wm1r7by/kPW8srF+/fkJff7LIoFRxyFEMMihVHHIUgwxKFYccxSCDUsUhB7BVuVjLfWNWXy5RSt1r47QgUGY+rmSc3G9Lp1fxxJ4LDEXilAW84/GSGo1GMyUpmphIFrYBa8zHq4CT4/GiixorUApOdw2Ox8tpNBrNlKVglohNfgFsEJHZwFuBm8bjRefVGcbNma5BmmZWjcdLajQazZSk6CwRZQbVzcd9QDOwCVinlOodj2vMrSsH4Gy3tkQ0Go1mLBS7JYJSqhvInW7lkIbKAGV+L2e6hzL+vWsgwn89fwSlFPetW8L06onL4tJoNJrJTNErkYlARJhbV8aZDDGRWDzBH3x7Mwda+/CI8PKxTn5x3xuoLJmSU6XRaDQ5KTp31qViVm0ZF3pDo47/dNtZ9l/o438+cC3f/6MbONoW5FsvHS+AhBqNRlP8TFklMr2qhPb+8Kjjj7x2hstnVfPWK2dyy5IG3nLFTB7aeEI3bNRoNJoMTGkl0hEMk0ikKvbPdA2y60wP77p6NiICwB/espD+cIxn9rUWSlSNRqMpWqa0EoklFF2DkeSxV491AnD78unJYzcuqmd2TSlP7b1wyWXUaDSaYmfqKhEz46qtL+XSeu1kF3XlfpZMT/Ww8XiE5uXT2Xikg0gsccnl1Gg0mmJm6iqRqhIA2vpTwfXtp7u5bkFd0pVlcduyRgYicb0jokaj0YxgyiqRxqQSMSyRaDzBqc7MFeyrF9QBhpLRaDQaTYopq0TqKgIA9A5GATjVOUgsobiscXQ75mmVJSycVs72U1qJaDQaTTpTtoKuMuDDI9A7ZCiRY+3GToeZlAjA1fNq2XS8K+vrJRKKLz97iCf3XGBd03TWVhbHPi0ajUYzkUxZS8TjEarL/Eklcrx9AIDFjRUZz18+q5rWvlDSchnJAxuO882WY1SV+vjOKyf56eFIxvM0Go3m9cSUVSIANWlK5ELvENWlPqpK/RnPbTL3Zs+0rW5fKMr/vHCUNy6fzq//Yg3vv34ez52KcapzYOKE12g0miKgaJSIiNSIyFMi8qyI/FxEAubxB0XkVRH59HhfM12JtPWFczZaXGYG3A+1jlYiT+9tpT8c477blyAi/M2dy1DAj7acHm+RNRqNpqgoGiUC/D7wVaXUm4BW4C0ichfgVUrdDCwWkaXjecFhSqQ/xIzqkqznzq4ppTzgTbq90nl89wUWTCvnmnm1AMyoLmVlg5efbz83rCJeo9FoXm8UjRJRSn1DKfWc+bQRaMPYS8RqA/8sqV0Ox4XqMj99phK52BdmelV2S0REmFdXzpkRe5BEYgm2nOhkXdP0YfUl18/00tYfZt/5vvEUWaPRaIqKgmVnicj9QFPaoReUUp8VkZuBOqXUJhH5U+Cc+fcu4Nosr3UPcA9AY2MjLS0ttmQY6A7T3hdj/fr1XOwdItTTlnNsWSLEwTMDw8452hMnFE1QMXiBlpb25PHFZWFAePiZzbzzskDW1+wcSvDcqRjTyoR183z4PJL1XKcEg0HbczGRFIMcxSBDschRDDIUixzFIEMxyeGGgikRpdS9I4+JSD3w38DvmoeCQJn5uJIslpNS6gHgAYCmpibV3NxsS4bNoYNsPHecq294A7FnnmP1FUtpXrMo6/ktfft4bNtZbrvttqTVcXTDceAAH3rb2mQBI0BLSwvLZ3pop4Tm5hszvl7XQIR3/PdGzvUY1lC4fDpfec8qW7LboaWlBbtzMZEUgxzFIEOxyFEMMhSLHMUgQzHJ4YaicWeZgfSfAn+vlDplHt5GyoW1Cjg5ntesKfMTS6iki6qhMrvFADC3roxgOEZPWprvodZ+GipLhikQi2vm17LrTA9KZY6LfP2Fo7T2hfjlfW/gvnWX8di2s2w63jmGd6TRaDSXlqJRIsAfY7ir/kFEWkTkfcAvgLtF5KvAe4EnxvOCFQEvYMRDAKpKcxtmc+sMo+hcT2pb3cNtQZpmZi5QXDW3lr5QjBMdo4PxwXCMH285zW+vms2qebV89PalNFaV8M2WY67ei0aj0RSColEiSqlvKqXqlFLN5s9PlFJ9GMH1TcA6pVTveF6zPGAoDWtzqsqSzDUiFpa10RE0zk8kFEcu9rN0+uh+WwBXza0BYP+F0cH1lkNtDEXjvO/6eQCU+r28d/VcNhxpz7hZlkaj0RQjRaNEsqGU6lZKPaqUGvddocpNS8Tq5FtR4s15fkOloUSsRb6tP8xgJM5l0zNbIpc1ViICR9uCo/723P6LTKsIsHphffLYO1fNIaHg2f253+q2U93c/pUW1vzbC6w/2JbzXI1Go5lIil6JTCRlSSViWSK53VmWEukIGi1NLLfW3NqyjOeX+r3MrSvLqEReO9HFTZdNw5uWjbVsRiWzakp55Wj2uEjXQIR7v7+VcCxBqd/LR36wjePto19fo9FoLgVTWomMdGdV5FEiFSU+ygPepDvrQq+hRGbVZq8vWdJYOUqJXOgd4nxviOvm1w07LiLcfNk0Xj3embVI8TuvnKRzIMK3PriaH/3pjXg9wlefO5xTbo1Go5koprgSMSyRdpuWCBjWSFKJ9BhusFk1mS0RgMWNlZzsHBiWobXzdA8A1y6oG3X+DQvr6RqIcLprcNTf4gnFjzaf5vam6ayYXc30qlL+4KYFPL23NSlTLk50DLD+UBuxuN6hUaPRjA9aiWAoEZ9HKPHln46GykBS6ZzvHaIi4KU6R1bX7NoyQtHEsLTgI6ZlsmzG6FjKlXOMYHymSvedZ7rpCIZ51zVzksfefd1cYgnF47vO55T76b2tvPlrL/Hhh1/jfQ9sIhSN5zxfo9Fo7DCllYgVE2kPhqko8Y3aFjcTDZUldJoxkQs9IWbVluUcN8d0daWnBR9tCzKntizpTktn6YxK/F5h7/nRiWjPH2jD5xFua2pMHls2o4pFDRW8eLh91PkW3QMRPvl/u1k+q4r/984r2Haqm/9+4Uje96rRaDT5mNJKJOA13n4klrDlyoLhTRs7gmEaK7M3bQTDEgE4P0KJLMmS0VXi87JkehUHMqQFbz3ZzZVzaqge0a5+7dIGNh3vIhzLbF38ZOsZegajfPGulfzhLQt5+8pZfPeVU8n3kYv2/jDrD7URDMfynqvRaKYeU1uJpLmv8qX3WqQrkZ6hKLXluWtLrHjJhV4jfqKU4nhHMOsOimBsjHVyRIFiLJ5gz7lerjY7Badz8+JpDEXj7M/S7PHR185ww6J6VsyuBuBP1y4mGI7x9N4LOWXfe66XO7/2Ih9++DXe9NUXaTXfg0aj0VhoJWLixBIZisaJxBL0DkWpKcutRKZVBAj4PElLpHswSiiaSFa/Z2LhtHLOdg8RTQuAH74YZCga55r5taPOX2kqlr3nRrvATnQMcLxjgLddNSt1/twa5tWX8cSe7PUo8YTi44/tpsTn4cvvXknPUJRP/XxPzveq0WimHlNbiXjTLRGbSsS0PPpCUVtKxOMRGitLaDezp6y7+Vk12dOCF06rIJZQnOtOucCOmDsqXj6retT5s2tKqa8IsPvsaCXy4iGjGHFd0/TkMRHhTStmsul4Z9YA+0uH2zlwoY9P/dblvGf1PO5bt4QXDrZlVFQjeWTLaX73m6/wb08f1JlgGs3rnCmtREQEv9cIipf67bmzrHjExb4QkVgiqVRyUV8RoHvACMa39hmKYUYuJdJg7PN+Im173WPtA3gEFkwrz/g+rphdnTGja9vpHmbVlDJ/xLibFk8jEktkVDwAP3ntDA2VAd56pWHB3H3zAgI+Dz/bfjbXW+Xx3ef55P/tobU3xDdbjvHFpw7mPF+j0UxuprQSgZQ1ErCR3gspi8WqEclniQDUlvvpMlN8W3sNi2Rmjq1459UZC356MP54e5C5deWU+DIruyXTR9ejAOw+28NKs4dXOtcvNGpUtpwYXR0fiSXYcKSdt1w5Mzkv1aV+mpc18sTuC1m7EkfjCT7/xAGumlPDix9v5vdumM/Dr5zk4kB+a+Rgax+f/NluvtFydJgbT6PRFDdaiZiLZLprKxdW51+rWt2OEqmvCNAzaFkiITxCxtbxFg2VAUSMfd8tjrcPsLixIuuYxQ0VDEbiyRYuAANRxanOQVbOrR11fm15gEUNFew9N9p62X66m4FInFuXNg47fseKGbT1h5N1LiNZf7CNC70h/uqNS/F5PfzNnUvxCLxwJncW2OnOQd77v6/ys+1n+dLTh/jHX+zNeb5GoykeikqJiEi9iNwpIg2X6pqWErHcWvkoNy2Rc6YlUluWew8SgLryAF2mO+tib4iGyhL8OZSWz+thWkVJsjEkwJnuQRbUj3ZlWSxqMLK90veAPxc07uhXZIijACyfWcWhi/2jjr92ogsRuOmyacOO37zYeJ5tz5Nf7jxPQ2UJzWYdy/SqUpqbprP5Qpx4jr3m/+2Zg8QSihc+1sy9ty3mkdfOsP10d9bzLQ5c6OPuBzfzZz/YNqwOR6PRXDqKRomISB3wOHADsF5EGs3jD4rIqyLy6Ym4bkqJOLNELFeTXUukPxQjGk/QNRihviK/4plRXZLc52QwEqM/FMsTRzEUTPreJZYbaVFDZgtm+cxqTnYOMBgZXgOy51wvixoqRtWjzK0rY3ZNKVtOdI16rXhCsfFoB+uaGvGlzeXbrppFT1hlTT9u7Q3x1J4LfPDmhcyrL+cvb19KVamPhzaeyPpewajR+f1vb2bf+T5eOtzOhx/eYrsKvy8U1e32NZpxomiUCLAS+Ful1OeBZ4BrReQuwKuUuhlYLCJLx/uilhvLthIxLREry8qOEqkzg+/dgxH6hqJU2xgzo7qUi33GNSy31oyq7EpkVk0ZHkm52QBaBxQ+j2RNJ26aWYVSo1vV7z3Xy5WzR8dRRIQr5tRk3B9l77leeoeirFk63Ii8xbRmXjnWkVGGX+06R0KR3FelosTH7147l2f3X2QgR4Hj1184Sn8oyiP33MT//P61HL4Y5IebT2c93+Lx3edZ/bnfcOMXfsPXddW+RjNmCrbH+kiUUi8CiMitGNbIZ4HPAY+apzyLsVXuqP98EbkHuAegsbHR0Yb3kZCx6LaeP0dLS/bWIRYDUcMtc7rdcLds37qJ4yWjFVAwGEzKceGCsRg+1/IK5zvCNJRJXhljwTBnOuK0tLRwqMu4w249cYiW/qNZx1QHhJ2HTtISMIoIz/VFaCj1sHHDSxnPb+83LJWnN2yla5bxVQhGFOd7Q5SFOjLKWB6OcKI9yjPPr6ckzQX43Ckj7hFvPURLy/CPaGaZ4onXjtCkzox6vZ9vGWJelYdTe1/D2hN5hlmH882ft7B65uivaCimeGTzIKtneDl/YBsAy+o8fPP5AyyOnszYhiYYDPLYUy/wqY1DzK30UFPi5SvPHka6TnNFQ+7MvGhC8cujUVoHErx5oZ+ldfYy+TKR/r0oFMUgQ7HIUQwyFJMcbiiYEhGR+4GmtEMvAP8CvA/oBqJABXDO/HsXxva5o1BKPQA8ANDU1KScbHhfv3cjZ/p7uWzRfJqbl+c9PxJLwPNPEcEPRGheu4ba8tHuqZaWFiw5/Ec7+MauzSy5YhWJfbtYPHcazc2rcl5ne+QQG84dZe2tt9G35wJs2cGb1t7A0hmZd1EEWLB3I5T5aW6+EYDPvPwUl8+bRnPzDRnPHwjH+MeXn6Fy1kKam5cAsON0N7zwCnfeuJLmFTNGjQk1tPLLY9uYseyaYdXzT/x0F9Mq2vidN68btYhftvsZDvV5Gfm5hKJxjv3mWe6+aQHNzSuSx9fEE3xj93P0lM6gufmqUTI8tecCofh2/urt13OTGafprDrLx366i/ol17AqQ1V/S0sLr/U2AGf4/p8101AZYN2XW2jpKOO+d9+ScX4s/uqRHTx+/DwVAS/7uqL84r4baJqZ/XOw2HO2l22nunjTFTOT7W/SvxeFohhkKBY5ikGGYpLDDQVzZyml7k3bCrdZKfVZZXAfsBt4JxAELF9M5YTIay54dt1Zfq8gQrKXlJ3UYKsaPhiO0ReKUl2WX3fXlAdQCoKhGG2mW2t6jrRgGO4CA+gJJ5iZI45SUeKjviLAma6UC+xUp9GC3oqxjMTqPHxshAts3/k+VsyuzmgFLKj20BEMJ9+Hxa4zPURiiaQisPB5PVw9v5atJzMH158/2EZNmZ/Vaa3077h8Bj6P8My+zFX4kbji/7af47eumsmc2jJKfF7+aM0iXjvZnXHTMIvXTnbxy53n+cs3LmX9x5sp9Xv416cOZD3f4jf7L/Kub7zMP/96P2/9zw22Nw7rGYzw9N4Lwz5HjaaYKZqYiIh8QkQ+aD6tBXqAbRguLIBVwMnxvq61+ZNdJSJitIyPxBK2x1lxlP6QESAfGbDOhNVevi8Upa0/TInPk7PlPMDMmtJkrCYWT9AfgcYccRSAeXVlnO1O7V1ysnMAEZhbl1mJzKkrQ4Rh+53EE4ojbf3J3lwjWVBtzNHIYkgrtrIqQx3LdQvqOHSxn77Q6PTgV452sGZJw7AAfk25n1XzatmcIegPcLArTn84xm+ntdF/+8rZiMCTe7L3EPv+q6eoLvXx582XMb2qlD9Zu5iWQ+0cbRud1WYxGInxD7/Yw/KZVTz2kZsRgX/61b6s51uc7hzkTV97iY/8YDtv+tpL7MvQyTnb9Z7ee4FjeodLTQEoGiWC4Y66W0ReArwYMZBfmMe+CrwXeGK8LxozlYjdOhFgWMGfz5M/NdiyRKwmjHYC69Y5vUNRegeNRo/5WtXPqC6lLxRjKBKnIxhBAdNz1KMAzK0v50yaQjjVOcjsmrKsFfwlPi+zqkuHjTnfM0Q0rlg0LXMW2JxKY25H3vHvP99HQ2UgY83Mqrm1KAUHLwxfrNv6QpzvDWXsIbZ6YR27z/ZkzNLa3R6n1O9JpimDoXSvnlebtY3+YCTG0/taedc1c5Lz8e7r5iICv96VXfH8aud5LvaF+czbV7B6YT1/sW4JG450sPtsT9YxAJ/+5V6GonG+/oFrKPF5+Niju7LucGnROxjld/7nFT7yg+285T9e4rn9F3Oeb9EXivLdV07y5J7sxaMajR2KRokopbqVUncqpW5VSv256drqA5qBTcA6pZS9WzMHpCwRe3UiQHLzqoDPY2sPEqtDsJMCRcta6QtF6Q9HqbJhvTQm94APJ90h+ZTI9KqS5J7xYOx7MidHc0iA+dPKOZWmRCyrZGRrFYvKgFBT5udU1/DOxPsv9HH5rMwusKWm2+zwiDqWnWd6ADIrkQX1ROOKPRn6ex3qTnD9wvpRyvHGRdPYfbaHochoxbP5eBeRWII702JDM6pLuX5hPc/mWKx/uu0sS6dXcsOiegDee/08Aj4PP99xLuuY3Wd7eOlwOx+9fQlvXzmbf3jb5Rxs7aflcFvWMQBffe4QR9uD/Pt7VrF8ZjUff2xXRustnXBM8f77N/FPv9rHn/9wO194Mr97DoyMxP/vsV184rHdw7op5EMplbTcNa8/ikaJZMNULo8qpbK3nB0DsYTplrLZ9gTcVLkblsh5s0Axn1sKoMpyZw3F6BuKJZ/nwrJeLBcY5I+jNFSWEAzHknfvHf3hvIpnfn15MnYCqTjKgiyWCBididPHKKU40TGQdV+VObVlVAS8HBmhRPad78MjcEWGFGTLnXawdfiYoUicc8EEqzJU7t+4yFA8O86Mjr+8dKSdUr+H6xfWDzu+ZkkDB1v7kl0I0ukeiLD9dLfpKjOUo9Uy5um9rVnv+v9v+zkCPg/vu34+AL911SymV5XwyJbRGW0WncEwP9pymveunsfvXjeXL/zOVfQMRvlxnlTnp05G2X+hj/vvvo7fv3E+39pwIq+VFIrG+b1vbeLnO87xfzvO8sGHtoyqL8rEjtPd3PyvL3DVPz+Tt/YnnVeOdvDAS8eGFdzaQVtVl56iVyITTdxhTASGWyJ28HiE8oA3efdmx51Vk6YQ+kNRe3EUM2DfOxRN/vPNqM6tEBoqjcwya4/2tv5wzpYsADNryugcCCc79J7qGiDg9eTsBzZ/WsUwJdI5EGEwEs9ahS8iLJlRxeGLw11gx8weYpncbbNrSqkq8XF4hBLZf6GXhCJjDzErk2tfhvYvO8/0sGpubQbrpR6lyFh0+fKxDpSCtcuG18vcuqyRC70hLg6OXuSUUjyzr5Xbm6YnP3e/18Obr5jJhiMdWYsoH999gWhc8aFbFgJw1dwarl9Yx0+3ZW+SGYrGefZklDetmMGbr5jJJ9+6nKpSH/e/eDzrGIAHN57gRMcAD3/oBh7+0A0cbQvyvVdP5RzTH4pyz/e34fcJ186v47OP7+eVo5nrhdL54eZTfODbm/nCkwd5x39vtKVIBiMx7vneVpo+/TSfeGy3bcvneE+cf/7VPp4/YM8NCIb3Ytupbq3gTLQSMT9YO7ENCysm4sQFVh7wJhdqO3uXWErDCsbbsUSSimcomuwanK86flqFoTA6gxEGIzGC4VheJdJYVYJSJFu5nO8JMbu2FG+OOVw4rZxzPak9UvK5wMDoB5YewIfcPcREhGUZWrnsN+Mq1v716dRXBJheVTLKeonGE+w/35dV8Xg9krED8ubjXVSW+Fg54lo3m0WXB7tGK4TjHQNc6A2NUjx3rJjBUDTOq1nazDy19wJNM6qGpRv/1lWzONoWzJoN1nKoncEY/MFNCwCoKvXz7uvm8tyBi/RncYPF4gl+sOkUa5c2sMb8Wbu0gYc2nsjZzua7r5ykvT/M13/vWh7+8PXMqS3jK88eyno+GDGvzz9xgLVLG3jsIzfTPRjlc4/nd7d95pf7+M2Bi9zW1MhPtp6xVUi67VQ3X9gS4juvnOSPv7uVx3IoX4tILMEHH9rC737zFW77UgsbjuSvLQP42bazXPMvz7HuKy1Jl2w+TnYM8OGHt3D3g5sdJVl8/on93Pv9rVm7RIwkkVD8z/rs9Wf5mPJKxCLXAjiSEr+zKncwFE+PuSOinbbzlUl3VpS+kL2YSE1aML4/FMPvIWvXX4uGqlQcpaPfUAr5tvxtNK0Xy2XW3h/Kq3hm1ZQRT6jk/vRWYH5eliwwMFxarX2hpMWTSBgusMUN2XeFXDq9clT68Yn2AQLe7Hu4NM2s4tDF4f9wRy4GCccSGRVPqd/LZY0VGbcw3nu+lytmVw/LHANDIU6rCHCsZ/Qd8qvHDCXxhsuGK5HVC+rwCGw/NdrVForG2X6qh1tHKB4rftNyKPPi9sy+Vqr8qU4CAG+9chaRWIIXDmaOv2w52cWF3hDvN11tAO9dPY+2/jBbT2bOhkskFI+8doY3LJnGqnmGNffHaxax/XTPqDhXOj/YdIqhaJx/+e0rWb2wng/fspAn9lzIGYM5fLGfn20/y5/euphvfXA171w1mwc2HE/eSGWT7zO/3EtNQNjyqTdy46J6/t+v9mVVpBbfbDnGxqMdfOzOZcyrL+OvHtmZd8yes718/LFdLKgvJxyN86ff25p3a+qBcIwPfGsTW092s+dcL3/40GvJm9BsKKX4ix/t4FsbTvDy0U7ee/+rtnrKffnZQ3z5mdzKPRd5V0ERmW/zJ3N+5yTBToDcwqk7CwzFY5nYpf7847weoarER18oSl8oZqu2JBkTGYrRH45hYwjTTEulMxihPWiY53YsEUi5wDqCERryKR5zjNWz6rTp2sqWSgzG/vTxhOKiOeZif4ihaJxFOboZz6svN11lKX/9iY4gM8uzJ0Esn2m4zdIzoXJtAmaMqR6lRGLxBAcu9GVUPCLCitnVnO4frUT2nO2lviIwaq+YihIfl8+qztiMcvupbiLxBLeMUDxz68qZXVPKtgxjlFJsOt7J5dO8w5TcdQvqqC71JZXZSF481I7fK9zWlOrqfPvy6ZT4PDyzL7MbaM+5Xs52D/G7185NHnvHqtl4BB7fdT7jmERC8bPt51izpCG5p84f3LSAeELlTEp49LUz+DzCvbdeBsB965YQiib46bbs8aRNxzvZd76Pdy3xM726lH942+X0h2P85LXsYwbCMb698ThvWjGDj75xKV969yq6BiL8eEvuGNSXnz1EfUWA7/3xjfzv3dfR3h/mu6+czDnm2xtOcL43xMMfvp5H7rmJroEw9794LOeY9YfaeOFgG59+2+U89VdrSSiVdz+fEx0DPPDScd5z3dyc5+XCzir4XeA75u9sP98B3uVaigJiuSm9jpSIcXfvNi3Y9gZYZX46ghEisYStmEhlwIdHUpZImS//e7IW//ZgOLnAT89TW9KQzAKLmL/DtpWI5Ue+2B+ittxPWSD7XFhZYtZdqJUiPac2u3xzR4wB4x9lZkX2uVgwrYJILDGsjf6JDqNeZn6WmM3ls6o53xsadkd5omOAUDTBFVnqZS6fVc25/sSo/VL2XehlRZYstWvm17L7TO8of/ou05V27fy6UWOuW1jPtpPdo8ac7hrkQm+I5fXD59zrEVYvrOe1LFbFS0c6uH5h/TA3bEWJj+sW1LE5w340YOyMCXDbspTiaawq4boFdbx4JHNc5GBrP+d6hnjHqtnJY/Pqy1k5tyZrzEIpxa92nae5aXrSdds0s4qr5tTw9N7suTg/33GOyhIfN5rtflbOrWXVvFp+lUXBAfzmwEX6QzH+ZO1iAK6eV8sNC+v56dbsbrDTnYO8dLidP7hpATVlflbOreW2ZY38aPPprOnbsXiCH285TXNTI6sX1rN8ZjXvWDWbH285k7PJ6HdfOcXsmlL+8BajmekHbpjPU3suJGvHMvGDTafwCHz8LU1Zz8lH3lVQKbVOKXW7+Tvbz+1Kqe+5lqKAJJWIA8eeFQtx4s5Ktz7sKpGKEi+tZlqwnYwuj0eoLvPTOxQlGIraUiJlAS+lfg+9Q1G6zY2z6ipyK6yk4ukPE4kl6BmM5lUi00dYIp02rBdLWVjbBF/stZIF8iuRM+aYaDzBme4hZpRn/6xSY9KKLjsGctbLLDIr+tPrZY6bHZSzZZxdPquKmBreaTkaT3C4NZhV8TTNrKY/HKN1RAX7/gt9zKkty7iz5tXzamntCw1L3YZUenSm3l/XL6znWPtAMs5lMRiJcai1j9UjMtQAblhUz/4LfRlTil8+1sEVs6uZNuIzvnnxNPae683oAnrZDLqP3Mfm9uXT2XGmJ2M23OGLQdr6w9x5+fAWPXeumMGOMz0ZXUCJhOL5g23ccfl0AmlxzTdfMYPdZ3uzus4e332BWTWlwzolvG3lLI60BbMWnz69z6gnenfanf7vXDOH1r5QxoxAMBI2WvtCvG/1vOSxd183l2A4RsuhzC7H9v4wG460c9e1c5Pr0u/dOJ9YQvHU3sw1TfGE4pc7z/PG5TPy3jjmwnFMREQqRMR9B7oixePAErHOdeTOSju31Oa4Ur83eXdsJyYCRlykdyhK0KY7C4wgft9QlGDIcAHlC/xXlPiSiQKdA4Z8DVV5AvhmHGW4Esk9ZlaNaVWYitRaSHNlgc2pNRb3s6YS6QiGiScU08qyf77z6q0xKYVwonMwawt9SLnhTg8r1DSUw4L6zOOsPV9GpkdH4omsvbiWmgrpyIgstQNmjU0mmsz+aiPTow+29uP3CrMyWGVXmS64g63DXXR7z/WRUJm7CqxeYGSp7T4zPOibSCj2nO3lugWjraQbF08jnlBsP90z6m+bT3SyuLFiVKueG8xsuEwBaUvx3LJkeOucNyxpQCkyts45dLGfroEIazIoK4BXMrj14gnFpmOdNDc14kmLn95hxqBeOpzZumo51E7TjKphbts3Xj4dr0dYfzBz3OrFw4b78NY0K+7mxdOoKfNnjVttONJOQsFbrpyZPHZZYyWLGyt4/kDmMXvO9dIRDPPWq2Zm/Ltd7MREPCLyARF5QkTagENAq4jsE5EvT0R79kLgRIlYQXgn2VnWHa3PI6OCrrnGWIuunYwuMBVCyL47C4yalP6QEUeBVF1LLqx9461gfD6rosTnpbbcT3syjhIedZc6kooSH6V+T9r+9CH8XsmZcTa9qgS/V1LWi9lGv7Yk+1zMMZsjpvcQO905kDNzzPpbuiVysnOQunJ/RusAUq6x08MKNU3Fk+VaSSWSliwQjsU53h7k8lmZFc+yLIWaBy70cVljZcZMxGUzjTGHRmSpWfUjmXbHXG5ef6TiOd4xwEAknjE2ZB3LlDm051wvV2e4zsq5tXgEdmRQPNtOdzOntmxUbO3KOdUEvJ6M8aTNZrbbzSM2XVs2vYqqUh/bTo126x240Ed/ODaqz9uc2jLm1JZldAWGY3G2nuxm7YjtEapK/Vw5pyZjijjAhiMdrF5Qn2yXBEY/uRsW1bPpeOYxG490UF8RGLUB3e1N09lyootwbLQbzFLAb1gytj0A7axm64HLgL8HZiql5iqlGoG1GJXkXxSRPxiTFEWAk+ws69xAnsyndCxLxK4ryzp30KykzhU7SCepEBwokeoyQ/EMhGNUBLzD7rSyjjGVVdISyWNVgJH1Ze2N0hEM02Bjc6768gBdA4bro60vzPSq0pxJEB6PMK2iJOnGsCr3cymRUr+X6VUlSUskFI3TPRhldo7mldWlfmrK/MNcYKc6B3IWXNaV+ynzGQrKwkowmJ/FeplWWUJduX9YX6xz3UMkFCzMcq3GqhJqyvwcHpGldqi1n+VZLJ5G8zojFc/RtmDW1jQNlSU0VJaMUjx7zY4BV2VQIjVlfubUlo1KSmjrD3GxL8wVGcZUlvhYNqOKXRkKIvef78voCizxeblyTnXGzLY95/qYXlWSvHmw8HiE6xbUsS3DGCsLbWThKRiW0msZYlAHLvQTiSe4NpNFtqienWdGt+gZjMQ42NrH9YtGX+emxdPMuNZod9vWU93cuKh+1P/udQvqiJjp6iPZcqKLphlVeW8A82FHidyhlPoXpdRupVQyIqiU6lJK/Uwp9bvAT8YkRQGxPng7C6dFUom4sETsZGYlx6S5vUpsusDK/F5C0Tj9oahtd1ZVqZ++UIxgKJZMLc5HdZmPvqFYsptxZYmdzbkC9AwZiQJ9oZitL29dRYBua3/63lDOrsQWlpUEJDsH15bm/qymV5ckrT67CQbz6suGWS9nu4eSrrFMiAiNZZ7hLrCuQcoD3pxKeE5d2TA/fb4aGxFhcWNF0r0GhmK80BtKutQyjVmaobjzeMdATrfe8plVo2psjrYF8YjhTsnE5bOqRikRa5HLFhtaNqNqlEsvGI5xsnMgY/cCMDoYHLrYP2px35/DFXjVnBqOtQ+MunM/2NpPfUUgY5r4yrk1dKQlplikrLgMtUZza4nEExn7ySUUo+qMgOTWC3tHFMb2DkY53TXIVRmuc42ZeDHSFaiUYt/53oyyOcVOYD13QrPNc4odBzokmcnlpso9X91GOunWh10LpjTgZSgSN2MiTtxZUYKR2DATOhc1adYLpPqD5cII+seSwdt87iwwFIJ1frfNrYWnVQbotPa07wvj9QjVgdxzMa2iJHkdKw7VmKfaf0ZVaXLhUErR1hdmRr706HJJxmvAsETm15fntK5m1wxXIpYLLVvmGBhulnNp17HqBebVZ++LNr++fFhcCIwkgFxK5LLGCk6mKSuAE50DzK0rzxozXDqjihMdA8MKFY+3G6+xNEtSwpLplZzrGRqWun2otR+lyNo9eun0KvpDsWFZd5FYgqM5Ok4vnVFFPKGS8lgcMK24TJ+TFc8aWeS652wv0yoCoyye9DEjFbDV9y2TQkiOGaGArULETLuRzqwpZXpVyah+cm39YTqCkaxK2wm2V0ER+U8ReVZEnhGRL4nI1WO+ehFgfY2dpPh6ku4s+0rEioP4nFgvaQrHiSXSORAhobBtiZT7TcUTilHlIPbSOxRlIGzcsdmJ2dSU+ekdjCRdYPV5ssDAsF4sS8RuG/10xXOxL0RjZUnemNe0ikAym6m9317zSkNZGe9lIBJnKBrPW2NTUyLJuBAYacuzMywy6cyuLeNCTyo760z3ECU+T86i0Dm1ZZzvDSXTSC3FlasuZ25dGW394eRdeL+5F30268V6vf5QjN7B1H3kyY6BZJ1HJubXlxNLKLpCKSVyqnOAKnN/m0xYGW/pi7tlaWVTcpmSEk51DhCNq2TywUgyxZPiCcXh1v68iQwj3XrH2oMsmV6ZUfEsnGYo2UMj4kmHLxoWT6YMxMoSHwumlY9SPNaWCtkUwtIZowtwLcWTyX3oFCfZWQeALwP/CbQBPxCRvxizBCMQkRkisiPt+YMi8qqIfHq8r5WOI3eWC0vECmY6aa+S7vqya8GU+b3J2oVSm5ZIecCIvQTD9t1ZVaV++kOxpCVSbiMYb2WOpRSPPYVguaaM/entBf0tJdI5EElmhuUiXSG0OaiX6QxGSCRU0iLJq0QCQs9gNFl4atTY5JZvdm0p/eaGZmBkkc2pK8v5nZ1TV0YklqDDfE/JDgE5LJE5tWUoRVJhWa66bEH/9NezYkNKKU52DLAoV1KCaUF1DKWUyMnOQRY0ZLfILNdYemzoVOcgItnf05JkUkJq0bVcgdne06KGCrweGaZ4WvuMItds7rlplSU0VAZGxZNOdGRv0ePzelg6vZJDI1x0+Sy/phmj2/qc6Bigttyf1bK/rLGSY+0Dw9x6+Sw/J9heBZVS/6uUek4p9aRS6ivAauDeMUswmq9g7mYoIncBXqXUzcDiicwEcxJYd2OJeJNKxEFMZJg7y9648rQxNjxMAJQFfAxF4mZg3Z4SMRRPjGAkRsDrsTUXNWV+BiLxpJIrtyFgXXmAvlCMcMzYVMqOJTKtIkAwbIzpG7LXvLK+ooRQNMFgJMbFvhBejySr+bNep7KEWELRZ96xgw0lYgb4OwfCJBKKzoH89TLJVGfTJdXRH8nbmma2OcZyaZ3tHsLvFWbkUIyWlWK5vi7226nLGZ5S3T0YZSASZ36OBAOr1U3bUKroMl9SQqqINGWRne4aZFZ1adYbrMaqEkp8nmFuvdN5XIElPi+zakqHufXsKOD59eXD4mPdAxG6B6M5W/QsnFbB2a7R7sNsCRNgKL8zXYPDFMKpzsGcc7dkeiXBcCyZqWiNqS71Zdza2ylu6kQ+IiJfwbBI7HX4sv/atwMDgFVq2gw8aj5+ltQuh+OG9Vk4iomYs+akYt2yQBz16PI6t0TSYyd+m9cqD3iJxI2iQbuWSHmJl4SCrmDEVjwEUnETK3PKjsKyCh+tf2Q7HZDrzMW/e8BIdbZjvVjWSmcwQmcwQl15IK91mt4B2a4SsWIzHf0ReoaixBMq7xhrER+W2ZZnjOUis6r82/pDTK8qzfmerIXaWnStpIRcnaDnJZWI8flctFHLM8ts1tlhdjROJJSRlJDD1VZR4qO61DcsM+l012DONGwRYU5tWXIOrDEVAW/O2Nqc2rJhyspOn7e5deXD+lSdyONqM8aUcbZ7KOlyNG5gwslC1kzMqy8nHEsMC+Kf7BxgYc5GpqYrsCNl9ZzMo7SdYNNrPowngTuBu4B/dXthEbkfSK+1fwFYB/wOxo6GABWA1TSnC7g2y2vdA9wD0NjYSEtLi205QmHjw9i5Ywf9J+wthhfOG2Naz5+lpSVzwVAwGBwmx7mzhntlaCBoW74zp1MVuptf3WgrG+zCmdSYWCRk61rnThuWQXt/iN6Oi/bGnDLGHDp1Aa9K5BxjzcW5M8aYLbuNfj67t7/GhRyV5ABnzxvusidf3AzA+ZNHaYnlbkF+odUY88yLL9PeG6bRN0SwLppTxvNtxpjnXnqV42ei+PK8J4BznYZb7jcbt3DG7Il1ePc2Wg9m/5z88RAgvPDqazSUGu+97fQxWqLZ39O5oPHar2zdReK8j9aeARaVh3PK1xMyxmzasZfyzkMcOR0ikFC0tLSM+m5aROLGYrZ59wGmDxxj01Hju3Rgx2aOZlE+SilKvPDaviMsiZ9mT7sxj+eO7aelK3tTv9oAXOyP0NLSQl9EEUso+i+eoaUle6uSal+cPcfO0tJi1Dccax3kqgZvznkoVUMcOD2UPGfH4RD1JYoXX3wRGP1/CuANhznUFU8e33gkggBHd2/hZJZ5iPdFONcd5YX16/GIsPmCMQ8Xju2jpS1zF+KhjiiReIJfPrsef2yQnz3zEgCDF0/R0pK5V1i3Ob+/fP5lltZ5iSYU57qHuK4+lnUeWgeM78L6TTuInDFuwg6dG2RhtcfRWpkN20pERB4F/kkpdQB4UEQeBnYAj7u5sFJqmCtMRD4DfEMp1ZPmFw1iuraASrJYTkqpBzC216WpqUk1NzfbliPw8m8gHOb61ddlLKjKxIv9++D0SS5btIDm5uUZz2lpaSFdjq3hQ3DiKHW11TQ3v8HWdfapo3DU+Ee8Y12zrbjNSf8JOLwfgKryUuzMxdnSU/zk0F7iCpoWZ39P6bRtPcMPD+wmHqigwadobr4167nWXPTtOs/D+3ZQ0TALjp7mjbetyZttFd1/kQd2b6Vh/jLYtocbr7mK5ityV9gGjnbw9Z2babryaiLbtrJ04VwqK9tzzkXtmR7+Y/vLLGi6ktKuU8zwR/N+TrNa+/nSay8xd8nlDLb2IweP8vY7cn9O7U++AAwxa2GT0Rvs5c3cesM1owrfho3pD/MPG3/DrIVLuOWG+Qw+/RSrmhbT3JzduxuJJfjrlqdonLOQ5ual/PuejSyYFqC5+YZR3810Kl58mtoZc2luXsFz3Xuov9DKnbevyzkPjVteoKKunubmq2nfega27ebNt96U2z21dyODoSDNzc1GseILG7jl2itpXjkr65glJ7bQEQzT3LyWRELR/+xTrFyW+/v6RPsuXjqS+uy/uPMlls0po7n5emD0/ykY/6ubXzzGmrW34vN6+NXFncyq6eSOHPNwruwUT5zYy+XX3sSsmjKObjgOuw7wjjeuyeoyUofa+N7+15h/+dX0n9hN5awV8PJrvPGWa7luweg6EYC5bf18bdtLNC5cTvM1czjVOYB6toU111xOc1qblHSGInE+ueFpamcvorl5CfGEouvZp/jdGxba+l/PhxN31g+An4jIHhH5DoabaTz3vLwDuE9EWoCrReTbwDZSLqxVwMlxvN4w3LQ9cRJY97hwZ6VXxNsN/Kfv0Gi3oj49jmI3xdca09Yftj2mwhxjmeLlNgooR24tbMedZbnkrPYvdmIi1n4twXCM/lDUVpaa5WrrHozSH45RWeLL+zlZMZH2YDjp1svnzqott64TSSvuzD0m4PNQVeKjy8xs6wyGk3vH5KJuWGZb/l0uwcxsG5EenS8pYVplCf0Rw/Kx6wqcVVOWbCaYdAXmiyfVGhlnVtPLzoFI3utY3aOt99Lal78+yUrjtWJDbf1hSnyenFthz7N6tnVZ7sP8czd3hPswWdOUw31YFjC6RVj/Q10DEWIJlXVrBKfk/U8RETH3O/8V8CsRWYWxoHswXFvJc8YiiFIqeSsrIi1KqT8x28tvEJHZwFuBm8ZyjVw4USJudkO0YiKC/es4CcJb+NPG2K1rTI+j2E0ltuIZ7f3hrKmPo8aUpMZ4PWLrWlbqcGprYRtKxBxj+eerSn2Qvflp6hxIFl3m8ulbpDYOM/qO2UlzDniFgM9DXyhK94Ax73VZ2qRY+L0eqkp99AxGk/ux2Mk4qzMz25RSdNjNUqtI1di094dyBtWTYypLkt2Z2/pCVJX68nZYmFYRYEdYmWPsKZGZ1aXJrtbW4pk3NlRTilLGoj6rupQuW4kMxnu20q87g5GcGWpAUslY76W115i7XPU/jaay6AiGqSXV4TrXPJT6vVSV+JLp6CmlnX/urKy75NyNsVLdwlbbExH5qIjMB1BK7TI79v4IuFJEvgv84bhIY6KUajZ/92EE1zcB65RSvTmGubuW+duJhZAw9aXdBXfY6zsI4DupKbHwp6X12hUvPUHArmJMXyQq7QbWLcUTDFMe8Nraw8VSPClLJP9C7UbxVJnpxkGzZYwdhVDi8+D3Cv0ho3LfUX+zoRgDZksbO8kMVr2MZb3YaTNTV+6ny8yWisQSebPNYHhKdc9QNGkF5RvTmbao2bFeGqpK6IsolFLJupl84+rN99w7lJYNl2chTCVZpBIZ8ilTy8VqdQ020sTzyGaOsSy/i32hvFtTV5f68HslqRAu9oWpLffnLSxOL6Ztt6lEZqclGNi1gO1iZ8V4C8Z93I9F5LyI7BeR48AR4PeArymlvjMu0mRAKdWtlHpUKZU94jYOOGkFHxuTJWIfN5aIb5gl4twFZjdtOf08OzUikErpbe8P204lthZmq4OvnYXaWpSdKJ5SvwefRwyrIhyz1TVZRFIdkB3U2FSXGZuNBcMx/F6xlXlXV+6ne9DYoAygpsyeJdIzGKEzaBV32htjubPsFndai5pSyn5XgYoAcWVsoNbeb9xU5HOL1iZ37kxTpvlcgeaYnsFomgLOo3jMGEbXgFED1DWQv5YnOSZNmeaz4kTErGmy6pNC9tyHlSXJz7StP4TPI8nrZ2NGdcpatKuA7ZL3W6+UCgHfAL4hIn6gARhSSvWMiwQFxnLCOdnZMOFCiTixdCxcWSJe5+6s9NiJ3bTl9PMqbDaHtCy3UDRBea29MVbcxFrY7LR/KfN78UiqrqKq1E++vjwiQlWpLxlHsbOnvfHaPpeWiNnw0uaYWtMScdJmpr48wNG2YPKu1b47K4xSynZx57SKAJFYwownxWy5wJIbmw2E6bLparOsop5B+3U5SUtkMJJM4893LSvW1TMYpWcoSkLlV8B+r4fqUl+yu0J7f/6N2sCoT+oaiECDFYOyM3cBTnakYiINlSV5Y3G15QF6BqPDLL9LZomIyDIxV1ilVFQpdeH1okAMDIXgpO2JFRNxYr04qVQfy5h0hWBXiaQrBLuWSLorr8RmX6/0O267loilNHqHoojYcyGKCGV+b7K4ys7dNBiuM2uMfSXiT8ZE7I6pLktZL3bnwbBEIo46BFhWxYCTJpkVAULRRDL4asciswL2ncEI/TbnIb0ux0hkyH+d2jLLzRSlYyBMwOvJmwCRVDxD0WTwP98deGWJD59H6Eqz4pz0eYvFDYVqxxWYHoPqHAjb7K5QMqy7gh1lUFvmJ5ZQDEbidPSHKfPnt/zsYudV/g+YJyKHgT3Abuu3UipzkcQkxImlkFIiDiwRc6F2oKts7zuSTrolYlcJ+V0okYCLLLDhLjB7isfnETyC0QvMby+OAobysayX8hIvmTdxHU6Z35u8S7OrEIa1jLH5T1lZ4uV8j9EhwLb1UuYnGIoltwawY/1Vl/oZjMTpD1mKJ/8Ya6G2sozsKGArA8nYxyZqa+5qkq4poyDUjivQGtMzlFLa+b4PScVjLu6QXyGISDIpwVrg7WxbYCltq7O13T5vRssYse0+bDCvE08oegbtxa3S584oVh17pbqFHXfWlSJSAqzEyMYaAN4BXCEiKKXGti1WkeBkcY+rS2OJ2I1pDLtO2oJuV77AGGMidt166VaE3bsgEUnuq+Kojb7ZiNJ6bIeygJf2ZOzFnvVSFjCUlZGdZXOM38dgOMZAOG672r884GMgEmcgEqPE57F1gzGyQ4AdJTIynmRHIVjnWPvY2LFeqtIz28L2suFqku6siG1XYMDnobLEl9z6OV3eXFiWn5UwUG/TFXi+J0TfkKlEbKSjT6s0khKUCjiw4kpIKGMeguFYxi7BI6lNm7u+UCxn6rFTbP0nK6XCwGsiElRKfdQ6LiKjd1uZZKT2WHduibjZDdEJbsaku6bsLrnDLBEXMRG7SiR9TJmDzblKfB4GI3FHY4Y3r3Swk6S54Npu/xLwcrY7TjBiP7BeUeJlMBp3FHupCHiNfViGorYVcGVaSrVxXTuZbcYcJzPbHNTldATDpgvMmeKxm5RQVeLDI6n6HyfbFvQMRvD7jPRqO9/XuvIA3QNRgmFD+dixEOrKA+w9l9pz3pZCMPu8DUQDxG26D1N1Q04sv1Rmm910dLs49ZcMqwVRSmXeaX4S4aYVfCJpiTip+RhbkNz2ddLG2H1LbmIiw6wXm3J6PJI818nWwlYsxemukJke56LM7yVqtv6w624r8/voDIZRyn6qc5npAnPizipPUwh2ZbPGWLULdnbHtKypZGabrcXdGGMpHnuWiKVE7C9qHo+YCsFMfrA5d3UVfnqGogyG47aTQOorAnQNRpIdp+3MubWBmqVE7CgeS0l3mW1qnClg+3OXdB+aczeeSsROseHXMdqbbMdZhuqkwkkr+GRMxIUl4qTY0FVGV9oYu8PTa0tsWyIuYiJgWAWReMJRvMeyKtwqEfuWiIudJAOepKvErjur3O8jEkvQ68CqSK/2t58ebY6x3Fk25s9ygZ3vta8QRrnAbNXYePF5TBeYg/ToVJZa3FatDBhxkZ7BCLVlftvzbW3UNhS1lEj+cRUBH+FYgm5zO2cnNU2d5t4q9pRIKp40EInbmrv0BINLrkQwAunXAB8EqkRkP7AP2A/sV0pN2q1x03GyXKcC6xNbfe4mxTd9cbetRNxYIunuLAdFlwGfB8LOLLOUJeJc8Xg9Yr+AMm2RtTsP6YuLm27GbiwRO1sEw/CuAgGbcRTLqrjQ4764066LrsxnpG5HYgnbVkVliY+BsLGPTb4qcovygJeOYJiBiP3Fs6LEx2DYsBb9ZpeBvLKNrE9y0F3B2qDLjvvQmt9WB4o+PT3aSWq5HewE1h9Ify4iczGC7FcBb2cS768O4KZbixsl4sIz5c6dlSaTXencBMnTs2LcbBPsaIdHUyE42b/F2hWy1MGYdHePm9b79q0XY0xCOVA8ab3KLrO5kVB6mxm7bhyrIDSZ6mzDugr4jEJNq5jNzqIGUOJN7fJod3G3XIH9Du6mK0p8pvswbtsVWBHwMRAxlJXdWNxY4kld5gZdttyH5vwmrUUb81Dm9+L3Cj1DEUeWnx3suLPmZzi81/xJ/3uP2aZkUuKk2DBu7UHiQImEoobP00lWhBt3ln9YTMRm6q2LtOBs4/OeaykRB5aZtaA7UVbW4m63hiV9jHHNibNe0sfYvSO0ZIsllO2F0FK+XYP5N7EaeZ3k7pg2rb8yvzcZe7FriZR4U/GaSpuKp9zMhnMUTzI3UHNqiSSU0fLEtssxaZENIWJvcbeUdGcyJmLfErGsRTsKQUSoKPHRa+6qWWnTJWoHO6/03SzHrVt4MR9/B/jeOMh0SXHTNTKeMD5wJwuutWFNrr2nR+J304BxjLUlTlKd3VzTWmidxFFSisf5GCeWSLoytKtEhicY2Fvc08fYvctNH2M3JmIp30gsYV/xmNcJhmMEfB77dTkBb7IAzu5CHXBhiZQHvJzuMuplnCzuA6ZrKtfOjsPHpGJQtq2XtD5vlYH8HZ3Tr9PlICZixbYu9jusafJ7k5l6l9QSUUrl3kxgkvO1917Nf/zmsCML4dNvW8FnfrmX1Vl6/mfi/dfPY/+FPv7ststsj3HivrFwY72kD3GStmzhRCFYr++oUNPjfIylcJxYIukuNvuxIedj3MSg0pWanW2FAUrSrAjbxZ1eD16PEE8oRwq41O9JZraV2LReAh5o73W2EJb5fXT0O1c8Q1Gj6NJ+IoPpCgw67/PW3h+2/RlZMbXesKVE8q9DPq+HgNeTUggOXIGW0r6kMZHXO+uWT2fd8umOxlw+q5qffuQWR2OmVZbwPx+41tEYJ12CLdy4o5y48jLhJLBuKQQnisfNGEshOLtOuiVibxEYe7W/8zFOLRGw3yQTDGtkIBJ3lA1XNswVaLcNTuqzsZN+DIZCsJpQ2q7LMd97p4utnNv7wjTOtucKtOagZyhquy+V5S7sM/dWsRu7KvWnlIj97gq+1JhCKhER+QDwTozOvgL8Win147EKIiI+4Lj5A/BRpdQeEXkQWAE8oZT63FivM5lw8k9s4SROM17jncRELH3lJCaSskScZ8M5sV7SOwTYVT5uCjVLxljc6cp6sbk4gfG9c6pE3MSTAmmnldpUPG42ULMUVCSesD3GsmD7HWQyWe/bSWGsNW9DMWO83Ru6skCqN5yT7gpOCk/t4iJniNuUUu9XSv2+UuoDpHYeHCsrgR8rpZrNnz0ichfgVUrdDCwWkez7gb4OuVSWSDpuRju5phlOcpSd5UtaIm4Uj+0hIwo1L00bfdsKIW1RctNmptzBopFMSnDkznKuREpcuALTr+MunuSs4zTYV8DD5sB2U1JP8sbKreVnN/mhPOBNbmNxSWMiGSgRkbcBZ4C5pPZAHys3AW8XkXUYDR7vxdiQ6lHz789iKKwjIweKyD3APQCNjY3jsvn8WAkGg2OWIxxPhf3tvlYkbYwbGTZvepW6UmfKa/eunQydzv4PkC5Hf9BIMDh18jgtLWdtvX5nh7kPQttF2++n9bxxxzVoXtvOXJw6YWQkCfbn+3BrLPl4x9YtnC3PPXfBYJDju3Yknx86sI/yzkN5r2P5zAFOnzpBS8s5W/L5PBBLQE9HW/I95ZuLeNSY72ho0PY8DJnpvT4PvPjii7bGSCKKddtiZ+4Azp2OJB8f3LcHf9uBvGOOnU99RhdOnxw2d9nm4mh3aivM3s52W/PQE07tFh4K9tmeO79ARIEkYrbHxCOh5OPNr748TCFnI9iTGrNv1w76jruxIUbjRon8OXAXRp3IGeAv3FxYRO4HmtIOrQfuUEpdEJHvAb8FVADWJ94FZAwqmLUsDwA0NTWp5uZmNyKNKy0tLYxVjkRCwXNPAth+rWg8Ac89BUBlZaV9GZ5+AoBbbrnF1n4Q6WOuvfZarluQvY1a+lyU73gJ+vtpWrqU5jWLbF3m5607oPU88+bMorl5pa0xm4YOwslj1NXW0Nx8i63P46T/BBzaj8L+fMcPXISdWwG4dU3+uWtpaeGypmvh1Q0AXLNqJc1N+WNyfaEorH8WgOVLl9C8drEt+cpanqE/FGPunNk0N1+VlCHX+6vbtYHWgT4a6425s8OPz2xlT8dFygI+23P3wwPPAMYCf9uaW3LuE25x1HscjhiKY/W113DzZdPyjhnacwF2bwfgyhVNNN+QqlrINhcN53ph80YALlswl+bmK/Jep3co9RnNbJxGc/MNeccAVLz0LJHBKFUVpbbnrmH/y5zp7wHgjnXNtly9v2rbydaLxnL6hptuYHGjvXqjfDhWIkqpQeAHY72wUure9OciUmI2egTYCiwFgqQsnUrcud8mLW7iE05asWTCzXAnYlp9x5y4wKz35CSO4nMTRxljrzI3LWPcdAhwMnelfi/9oZijBAPLPWI3QA4pF5sTl0z6d9VdZpvzLQjsfkbpLiK7wfj0MU47TkPUdlwIUu4sn0dsf8dLXHzv7DDmVxKRT4yHIMD3RWSViHiBdwG7gG2kYi6rgJPjdK3XLWMNrDvp7ZW8pgPNYzllnMREXAXWzdd3Mh1OFtpMY9woBNvB+PSMLgcLQLJDgKMEA+cdAqzF3UkcJf1UV0rEZl1OujK0+71LH2NXmQa87uIbqcJY55+rk/l2kwRiBzfZWY+mPwWuBv5tHGT5LPAj8zV/pZT6jYhUAxtEZDbwVoy4iWYCcaODnCiR5NbCTqwKF+m61t26k642bvqbuWqJfwnbzFjzMNFK2/psnCxq6SJNqBXnwhJJH2N37kSEEp+HUDThyKqw5szJGMsCdqPonY7LK4uLMX1KqT+xnojIN8dDEKXUXowMrfRjfSLSDNwJfEkp1Tse19Jkx03NiJMhbtrouylQdOeaGlubGbvvyU267vBrOlcITlxg1jw4GeNPLmoOFsK0t2738xq2/bOLHTXdZLY5uWsv9XsJRROOPlfrXCeWSEppO3c5pl9zPHCjRD4/4vk/jIcg2TD3LHk074mvUx7+8PW2WzWMBxNtibhxZ6VSfF1YIrZHuLNE0v8x7bfed76oZbtmPpJKxIX14swS8Zi/nX9GTnDj13fTcXqYC8xRx2nTfehiHtwoBEfKKt31Wgh3logcwUi93SUiO4FdSqmTSqmucZNGM4p1NjJ3xhM3loiTtdeyRJwshCLOF7XkP74Td5ar1vupMXbnblhA2VWnZvtjLHeek62WLYvPmSVixaCcJEzYPjXtOmNzZ7mxXpwpYEuZOh/j5LtgjXHiPnT7nvLh5JXuB1qBToz4xF4R2SMinxWR8duwV1NQ3GVnuSg2nOCdHr3mucqBFnETWHdjvaRP11gVVz6SO3e6uDN24z50puiN306+c646BLiw/NLfhxv3oZskECefkfUdcBsTGU+cuLP+QCl1tfVERP4X+DDQB3wV+GiWcZpJhJsGjG5cYE4WT3WJtiN2oxDcNLxMn2J3DS+dWCLOExm8Y4iJONvt0/jtJPlhrIF1d3v0uEhkcJFa7mjMGDLoxhsnr9orIsnAt1JqJ3CTUuorwBvGWzBNYXCTIOxoL5aEVSfiYCE0fzv7J3OhRFxYBW5SqtPTqN1Yfm7cWW589E7ujP3JGJR9jeDAoEpdx0U23PAxE/y9ME91VdPk4MvgKvlhHIPpw2RxcO5HgB+Y8ZCdGNXmVp2/vc2ONUWPO0vEgTvLKjZ08I/ppkDRen0nd7lu7tTcJSKkHrtRIk7cWan5dhGMd7FAJRzMtxsrLt01ZVe+9NPcxKBcWS+u3IcO3FmmTE7+96zA+hjrkUdhe3aUUgeAG4CngenAUYxeVxXAI+MrlqZQTHTFesqqsP+PGYkZ9ypOCris13eWnTWxCtQi3XJzU9zpxopzElh3Ve1vnht3oEXcWCLpi7NdC9htjU3ymi6+4G4y21zFUVxcx8WO4LllsXuiiNQDf4OhQPYD3zPTbwGmVIv21zMTHVhXLiwRa2thZ0rEskScBNbHFiS3PSbtsRtLxsnCoVwsam589JZv3sl8u3HRjzVm56oWaIK7aY9lqwO3SQnjiZNXfQToB34NlAMbRcRehzHNpMHNnbGzYkPjt5NAbyhqdFR10o8ofQ8Ju7jbFXJs8+UmpdqREsF5SrWb7Czr9Z24s9xYIm48Menz5cqdNdHxOxdWhaW0E/a/3kWRnTVLKfUl8/HjIvITjDYluhXJ64gJb3viItMqZLqznGx1W21udzwYjuc5M4WboKubhpfD3VnOcTJ3bvZvcVXlbp6bcGSJTOwNi0X69/NSdTLwulLaEzvfbr7fdnAyo10jsrOOY1gkmtcRrooN3fTOcuTOMi0RBxW91eZe1cFwLM+ZKdyk+LqxRMY63pX70MF7s17fyeJp3Rk7iYm4uWEZa1seN0rfVWNOF+5DJ7KlLD8n8z0xSsSJJXIv8JiIbMCoXF8BHJsQqTQFY6zZRvlIpZzaX9TCLtxZ1WXGV3swYt8ScXPHKWP1EEx0TMT87cSqsM50U7fgJGjrSok4HzJs8XTzeTlLqXZuaVsKx1FKtdd5IsME6RBH2VkHMTaFWo8RXN8F/N7EiKUpFO7u9CY2XfePzQ2YnGyiY1kiAxH7loi7FN+xWiLOxzi5Y7W2Q3XS3C95HUcxEefuFTczN1ZL2Y0l4qgw1sUYa56d1YmYlp8Dpe1m7mzJYuPC38OoC9mJ0S9rQpshisg3gKeUUr82nz+IYfU8oZTSWWBFiKsUXwf/ZO9cNZt3rprtSKZSv5fLGiv48+Yltse4SfEd+yZgbvzt9sdYd6puAspOhrhzrziVyK0lkv7YjTtrYqvc/S4KB1O94SZ2vm3JYuOc72JsCPWHwCoRqQIOkFIqPx0vYURkLTAzTYHcBXiVUjeLyEMislQpNWqPdc1wPnbnMiMdNnH6klzPVbGhi/iDU57/WLOj890EXcd6c+cqsO7EEjGz05y0x7Be3kmmnpvsLHcxETdjUoPcfO2cLO5Jd62TrQ7Mc510P/AmA+u2hxQuJqKUeh543nouIj4My2AVcCMwLkrEbOL4LeBJEfltpdQvgWZSbeCfxdjlUCuRPHz0jUsBaGkpRiVi/J6oTJGx4EamggTWHSyEliXipM24haMWJi589O4skUs/324sEScxEetMZ33HLk1mmx1sBdZFZDnw28AcDI/EeYzdB7/v9sIicj9G6xSL9RhFjF8CPioi84EK4Jz59y6MmEym17oHuAegsbGRlpYWt2KNG8FgsOByuJHBjcwvv7yRcn/2b2i6HHHzzvjVV16mzE2hgEvszEUsbQG0Ow9xh2NGyrFhw0vD9nmww6ZXX6U6YG9MJGYkFuzcvpW2w56MMozkzJkwAIePHKMlccbWdY52G9cZHBqyPXehoSGsJdTumPbBVGGEm+/qhpdeHKZI7Hwvtm7ZxLFSe4pkaGgIgAP79xJoP2hrzOnTEQCOnzhBS8u5PGcbHGo1Yn19/f225+FgVyrJZDzXJjsxkU9gBNAfAbaYh+cCPxaRR5RSX3RzYaXUvSOu83XgAaVUq4j8AGPzqzNAmXlKJVkSAZRSDwAPADQ1Nanm5mY3Io0rLS0tFFoORzI8/QSAM5nNMbfeupbKkuxfpWFyPPskKMW62251VIE+VuzMRSKhDPmwPw9KORuTlMOcu9tuu9W+lWDN95o3UFtur12devZJQLH2lpuYV18+XIYsxKZf5OmTW2lauoTmWxfbus60s72weSMlJaW25+7cr18AjEXX7piz3YPw0npHY4Dk3K1rbh7m3so5F+aYtWveQENlia3LlG55AYaGuGbVKm5d1mhrzI7oYTh+hHkLFtLcvMzWmNDeVti5jcrKSpqb19oaU3a8E7ZsAhzOXR7sWCJ/DFyhlIqmHxSRrwL7AFdKJANHAesbuxo4BezAcGFtwnCfHRqna2nGESeWtZvsrEuFq468Yw2su3DPOJEzmZ3lICZyx4oZPHrvzaycW2N7jM9NdtYlqhMZ63gnFevxuPMOAVa8Kuqgu4LPTUxkgv7n7CiRBDAbY1FPZxapLr7jwYPAQyLyfsAPvBujzcoGEZmNsRGWro4vQpz4mf/xbSv4lyf2u/Lpvh5x1wre+SCnMZEbFtU7Ot9aNCc+JnLpcZJJ2Dtk3GtPq7Tf2NzKnLMajdrBm+xSPTmys/4aeN7cHtdykM4HljCOG1EppfqB94w8LiLNwJ3Al5RSveN1Pc344WQh/KM1i/ijNYsmTphJxkRXrFs4yc5yQ6pOxP6YS5WdNVacKJEBs7h1WoUDJeJzoUQmuJ7LCXays54WkWUYbeDnYNwMnAVeU0rZLwd2idkpeEJrUzRjY6y1ElMZVym+LlbfiVcizrv4uio2LIAt4sSdZVFnM2YF7pRIqku1fZkK3fZkGUam1DNKqaB1UETeopR6ekIk00waJurLORVwM3VulPZEuw/dxEQmiyXiJJbw5Xev5Fe7zjsaY7mznMRErM/TSRp2wdxZIvKXwH0YBYYPishfmTUcAF/A2KRKM4V5PemQ25dPZ82Shkt2PVdtPBysBo9/dA2vHut0fA2nWHfrE9+A0fmYS8l7Vs/jPavnORpj7VcSdhJYd7Fz50RZcXYskT8FrlNKBUVkIUYTxoVKqf+kMHEuTZExUb7WQvDQh64vtAjjypVzarhyjv0sK7d4LpF7pRDurImm1NWGXs537ixksaHXcmEppU6age7HRGQBWoloNBpSKcRrltq34l6Plogb1i2fzhvn+/jnd1xhe4zl0pwsreBbReRqpdROANMieTvwEHDVhEil0WgmFaV+Ly1/18zMmlLbYyZLiu9E4/d6uHtFCdOr7c+dq104Jyi3ws7LfhBoTT+glIoppT4I3DohUmk0mknHwoYKR10ILtWmVK9HkmnHkyE7Syl1NsffXh5fcTQazVTBzY2xViEGqS6+hc/Omvh+3BqNRpMBd21Pxl+OyYhlVTgJrE+UCnayPa5Go9GMG8XuzvrFfW+goz98ya7nBGsWJno7YjtoJaLRaApCsWdnXT2v9tJdzCGSDIkUPjtLu7M0Gk1BcNf2RAOpepnJ1PZEo9EUmIc/fD1HLvYXWoxxw41rSmdnGSQtEScV69qdpdFMbdY1TWdd0/RCi1FQtAoxKAsYqdSXz6qyPaaQ+4lcEkTkz4D3mU9rgc1KqXtF5EGMPd2fUEp9rlDyaUbz/Mdu42THQKHF0EwhtCFi0FBZwiP33OSopc1ETV3RxESUUt9USjUrpZqBDcC3ROQujLYrNwOLRWRpQYXUDOOyxkreePmMQouhmUK8HntnueWmxdNybks9kikTExGROcAMpdRWEfkgqb1EnsXYKvdIhjH3APcANDY2jusm9G4JBoMFl8ONDBMh82SdC6fYef2RchRiXorh87DksO6P7coTiaeCAG7ew8gxxTQXEy1HTyjVJXg8r1UwJSIi9wNNaYdeUEp9FqPt/DfNYxXAOfNxF3BtptdSSj0APADQ1NSkxnMTere0tLRQaDkcyfD0EwATIvOkmwunOJi7pBwTON+2ZSgwLS0tfPDmBtYubaR5hT2LNhyLw3PG7hOO3kOW+S6muZhoOdr7w9DyG2B8v3cFUyJKqXtHHhMRD7AO+AfzUBAoMx9XUkTuN41GM3Y++9tXOjpfu7PcM1XanqzFCKhbNus2DBcWwCrgZCGE0mhyUV8R4Nr5tYUWY0qgA+vuKdge65eYNwMvpT3/BbBBRGYDbwVuKoRQGk0utv/jnYUWYcqgdYh7pkTbE6XUp0Y87zM3wboT+JJSqrcQck0FXv372xmKxAsthkaTE11s6J6pYomMQinVTSpDSzNBzKopy3+SRlNgJupueiowJSwRjUajyYXbu+kn/3LtOEsy+ZgydSIajUYz3qyYXV1oEQqO7uKr0Wg0GtdMVDhJKxGNRqOZAmglotFoNBrXaHeWRqPRaFyjlYhGo9FoXDNV2p5oNBqNZgKYqGJDrUQ0Go1G4xqtRDQajUbjGq1ENBqNRuMarUQ0Go1mCvEX65aM6+sVTdsTEakDfghMB7ZZm1aJyIPACuAJpdTnCiiiRqPRTGpOfvFt4/6axWSJ3A38UCm1GqgSkdUichfgVUrdDCwWkaWFFVGj0Wg06RSTEukErhSRWmAecAZoJtUG/llSuxxqNBqNpgiQ1E60l/jCIvcDTWmH1pvPDwJzgfuA/wX+Sym1S0TeBFyrlPpihte6B7gHoLGx8bpHHy389iPBYJDKysopL0OxyFEMMqTL8aGnBwD4zlsqCiZDoXErx3jO3WSfi/Fk3bp120xPkDOUUkXxAzwEVJuP/xZDKfwncJN57C7gU/leZ9myZaoYWL9+faFFKAoZlCoOOYpBBqVSciz4xONqwSceL6gMhcatHOM5d5N9LsYTYKtysXYXkzurDrhKRLzAjYACtpFyYa0CThZGNI1Go9Fkomiys4B/BR4GFgCvAj/GiNlsEJHZwFuBmwonnkYzfvzLu64kGIoVWgyNZswUjRJRSm0Brhh5XESagTuBLymlei+xWBrNhHD3TQsKLYJGMy4UjRLJhlKqm1SGlkaj0WiKiGKKiWg0Go1mkqGViEaj0Whco5WIRqPRaFyjlYhGo9FoXKOViEaj0Whco5WIRqPRaFyjlYhGo9FoXKOViEaj0Whco5WIRqPRaFyjlYhGo9FoXKOViEaj0Whco5WIRqPRaFyjlYhGo9FoXFM0SkREFonIEyKyQUT+Pe34gyLyqoh8upDyaTQajWY0RaNEgH8D/kUptRaYKyLNInIX4FVK3QwsFpGlhRVRo9FoNOmIsbVu4RGRncANSqmIiPw38BvgjcDTSqknReT9QJlS6uEMY+/B2JOdxsbG6x59tPDbjwSDQSorK6e8DMUiRzHIUCxyFIMMY5HjTH+CMh80lI39Hniyz8V4sm7dum1KqdWOB7rZmH08foD7gZa0n88AnwfeARwBKoEHgVXm+W8CPpnvdZctWzYOW9aPnfXr1xdahKKQQanikKMYZFCqOOQoBhmUKg45ikEGpYpDDmCrcrGWF2xnQ6XUvSOPicga4OPAd5VSQREJAmXmnyspLvebRqPRTHmKbVHeCcwHvmo+3wasMR+vAk5eepE0Go1Gk41i22P948BXlVKD5vNfABtEZDbwVuCmQgmm0Wg0mtEUlRJRSv3TiOd9ItIM3Al8SSnVWwi5NBqNRpOZolIimVBKdQOFT7fSaDQazSiKLSai0Wg0mkmEViIajUajcY1WIhqNRqNxjVYiGo1Go3GNViIajUajcY1WIhqNRqNxjVYiGo1Go3GNViIajUajcY1WIhqNRqNxjVYiGo1Go3GNViIajUajcY1WIhqNRqNxTcGUiIjMEJENI449KCKvisincx3TaDQaTXFQECUiInXAd4GKtGN3AV6l1M3AYhFZmulYIeTVaDQaTWYKZYnEgfcBfWnHmkm1fH8WY0fDTMc0Go1GUyRckv1EROR+oCnt0AtKqc+KSPppFcA583EXcG2WY5le/x7gHvNpWET2jpPoY6EB6NAyAMUhRzHIAMUhRzHIAMUhRzHIAMUhR1P+U0ZzSZSIUupeG6cFgTLzcSWGlZTpWKbXfwB4AEBEtiqlVo9J4HGgGOQoBhmKRY5ikKFY5CgGGYpFjmKQoVjkEJGtbsYVU3bWNlLuqlXAySzHNBqNRlMkFNP2uL8ANojIbOCtwE2AynBMo9FoNEVCQS0RpVRz2uM+jED6JmCdUqo30zEbL/vA+EvqimKQoxhkgOKQoxhkgOKQoxhkgOKQoxhkgOKQw5UMopQab0E0Go1GM0UoppiIRqPRaCYZWoloLikiUi8id4pIQ6Fl0Wg0Y2fSKhE77VAuRcsUm3KMavFyKWUQkRoReUpEnhWRn4tIoEBy1AGPAzcA60Wk8VLLkHbeDBHZMd7XtyODiPhE5LSItJg/VxVCjrTzviEi7yiEDCLyZ2nzsNOsKSuEHHUi8qSIbC2gDItE5AkR2SAi/z4RMpjXybseOVk7J6USsdMO5VK0TLEpx6gWL5daBuD3ga8qpd4EtAJvKZAcK4G/VUp9HniGLMWjEyyDxVdI1SBdahlWAj9WSjWbP3sKJAcishaYqZT6dSFkUEp905oHYAPwrULIAdwN/NCs1agSkXGt2bApw78B/6KUWgvMFZHm8ZTBlCPveuR07ZyUSgR77VDsnHMp5MjU4uWSyqCU+oZS6jnzaSPQViA5XlRKbRKRWzGskVcvtQwAInI7MIChUMcbOzLcBLxdRLaYd3wTkWqfVw4R8WMs2idF5LcLIUOaLHOAGUopVwVv4yBHJ3CliNQC84AzBZBhGbDdfNwG1IyzDGBvPWrGwdo5WZXIyHYoM1yeM+FyKKX6bKYmT5gMFiJyM1CnlNpUKDlERDC+xN1A9FLLYLry/hH45Dhf27YMwGvAHUqpGwA/8FsFkuODwH7gS8ANIvLRAshgcR/wzXG+vhM5NgILgL8EDpjnXWoZHgP+yXQtvgV4fpxlsLseOVo7J6sSsdMOxVbLlEsgx0RjSwYRqQf+G/ijQsqhDO4DdgPvLIAMnwS+oZTqGedrO5Fht1Lqgvl4KzAR3antyHEN8IBSqhX4AbCuADIgIh7z2i3jfH0ncvwT8BGl1GeBg8CHL7UMSqnPAU8BfwJ8VykVHGcZ7OJoXZusSsROO5RL0TKlGNqy5JXBvPv+KfD3SqlTBZTjEyLyQfNpLdBzqWUA7gDuE5EW4GoR+XYBZPi+iKwSES/wLmDXOMtgV46jwGLz8WpgvL8bdv8/1gKb1cQVrdmRow64yvxMbsTolnGpZQDYCcwHvjrO13eCs3VNKTXpfoBqjH+8r2KYnquAz+U5p6YQcqSd21LAufgzDPdRi/nzvgLJUQc8B7wEfAOz2LUQn8dEfSY25+FKDEtsD/D5An4vqjBuLl7CiE/NKcTnAXwBuGsi5sHBXNwA7MO4C38OqCzQXPw/4O6Jmou067SYv1eMde2ctBXrZpbBncBLyjDHXZ1zKeSYaIpBhmKRQ8tQXHIUgwzFIkcxyGAXJ7JOWiWi0Wg0msIzWWMiGo1GoykCtBLRaDQajWu0EtFoNBqNa7QS0Wg0Go1rtBLRaC4xInK/iLyh0HJoNOOBViIazaXnRozdOjWaSY9WIhrNOCIi7xSRx0Yc+zMR+S/z8eXAYWCRiLSLyEmzBXqXiBwTkepCyK3RuEUrEY1mfPk88M8jjh3DqAwGeCvwtFLqKEbTv7uVUldjVLC/Syk1Ud2eNZoJQSsRjWacEJFVgEcptVdEFojIn5l/8pPqxfRm4Gnz8RXAXvPxcuDQJRNWoxkntBLRaMaPqzGa14HRMsLqzrsC2CUi5UCtUuq8iJQBpUqpbhGZB3QqpSKXXGKNZoxoJaLRjB8eoNLsBHsXxg55ZcCHgB9htDtfb567AqO5HcDlaY81mkmFViIazfjxJEZr9Z3A/2K4q7Zi7NmxHTMeYp6b7soaAq4VkeWXVFqNZhzQDRg1mkuEiGwHblRKjfeOjhpNwdBKRKPRaDSu0e4sjUaj0bhGKxGNRqPRuEYrEY1Go9G4RisRjUaj0bhGKxGNRqPRuEYrEY1Go9G4RisRjUaj0bhGKxGNRqPRuOb/B3+3kH0CGpnTAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#导入使用的库\n",
    "import numpy as np\n",
    "from scipy import signal,fft\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.ticker import MaxNLocator\n",
    "\n",
    "#计算矩形窗的wn值和幅度响应\n",
    "N = 51 #矩形窗长度\n",
    "wn = signal.windows.boxcar(N) #矩形窗的wn值\n",
    "N0 = 2048;N1 = int(N0/2)\n",
    "Ha = np.abs(fft.fft(wn,N0))+1e-10;Ha = Ha/np.max(Ha)\n",
    "Ar = 20*np.log10(Ha);freq = np.linspace(0,1,N1)\n",
    "\n",
    "#绘制矩形窗\n",
    "fig,ax = plt.subplots();ax.stem(wn,basefmt=\"\")\n",
    "ax.set_title('矩形窗');ax.set_xlabel('n');ax.set_ylim([0,1.5])\n",
    "plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签\n",
    "fig.savefig('./win_rec1.png',dpi=500)\n",
    "\n",
    "#绘制矩形窗的幅度响应\n",
    "fig,ax = plt.subplots();ax.plot(freq,Ar[:N1]);ax.grid()\n",
    "ax.set_title('矩形窗的幅度响应');ax.set_xlabel('k')\n",
    "ax.set_xlabel(r'$ \\omega / \\pi $')\n",
    "ax.set_ylabel(r'$ 20log_{10}| H ( \\omega ) | $')\n",
    "ax.set_xlim([0,1]);ax.set_ylim([-100,1])\n",
    "ax.xaxis.set_major_locator(MaxNLocator(11))\n",
    "ax.yaxis.set_major_locator(MaxNLocator(11))\n",
    "plt.rcParams['axes.unicode_minus'] = False \n",
    "fig.savefig('./win_rec2.png',dpi=500)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "93f1e3fa",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
